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Abstract. Vibrational properties of graphene nanoribbons are examined with density functional based 
tight-binding method and non-resonant bond polarization theory. We show that the recently discovered 
reconstructed zigzag edge can be identified from the emergence of high-energy vibrational mode due to 
strong triple bonds at the edges. This mode is visible also in the Raman spectrum. Total vibrational 
density of states of the reconstructed zigzag edge is observed to resemble the vibrational density of states 
of armchair, rather than zigzag, graphene nanoribbon. Edge-related vibrational states increase in energy 
which corroborates increased ridigity of the reconstructed zigzag edge. 

PACS. 61.46.-W Structure of nanoscale materials - 64.70.Nd Structural transitions in nanoscale materials 
- 63.22-m Phonons or vibrational states in low-dimensional structures and nanoscale materials - 78.30Na 
Infrared and Raman spectra, fuUerenes and related materials 



In past decades carbon nanomaterials have shown their 
rich properties in several applications T,'2]. Here graphene 
nanoribbons are not an exception, and their use in many 
applications, among which transistors in nanoelectronics, 
have been investigated considering different edge structures [3l 
21 . Low dimensionality and edges make ribbons particu- 
larly fascinating both for theory and applications. 

The precise edge structure in graphene nanoribbons af- 
fects many properties like chemical reactivity [5], electronic 
structure [6] and vibrations [T. Vibrational properties play 
a role in structural stability structure identification i^S] 
and ballistic transport through electron-phonon coupling[9] . 
In structure identification scanning tunneling microscopy [10| 
lllj can reach near atom resolution but analysis of the full 
edge structure and properties is often ambiguous. In this 
case Raman spectroscopy j 1 2 ( 1 1 3 \ \9\ is a valuable tool. 

Vibrational properties of graphene nanoribbons have 
been recently studied for acoustic and optical phonons, 
symmetries, and Raman activity as a function of ribbon 
width [l2l[8]. In this paper we consider recently reported 
self-passivating edge reconstruction of zigzag ribbon [T?^, 
and investigate how the reconstruction affects Raman spec- 
tra and vibrational density of states (VDOS). Edge- loca- 
lized contributions to Raman spectra and vibrational den- 
sity of states change. It turns out that high-energy modes 
due to triple-bond vibrations, as well as overall changes in 
the vibrational density of states and changes in rigidity of 
the edges, make the reconstruction identifiable and visible 
in Raman spectra. 





Fig. 1. Examined structures of zigzag, reconstructed zigzag 
(zigzag57) and armchair graphene nanoribbons. Selected edge 
atoms are drawn in gray color. Coordinate axis is representated 
at lower right corner, where the ribbon direction (=periodic 
direction) is z-direction. W =width in non-periodic and L = 
length in periodic direction. 
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We use density functional based tight-binding method 
to calculate forces [1 5; . to optimize structure and to 
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Fig. 2. Raman spectra of zigzag, reconstructed zigzag (zigzag57) and armchair graphene nanoribbons (different rows). Columns 
denote different incident and scattered hght polarization spectra (uppermost symbols xx, xz and zz). Edge- localized modes are 
assigned with symbol (L). Symbols (1) and (2) are refered to in the text. Raman spectra are drawn using Lorenz distribution 
with full width at half maximum of 5 cm~^. 



calculate the vibrational eigenmodes for the systems. For 
Raman spectra calculations we use non-resonant bond po- 
larization theory [TTlllSj. that has been used succesfully for 
ribbons JJJ and for non-identical carbon nanotubes 19,20, 
?IT\ . The slight consistent energy shift in the high-energy 
modes [19] does not affect the qualitative changes in Ra- 
man spectra or VDOS we want point out in this paper. 

Figure [T] shows the atom structures of zigzag, recon- 
structed zigzag (zigzag57) and armchair graphene nanorib- 
bons. Widths and numbers of atoms of the ribbons are 
chosen to be easily comparable; also other sizes were sys- 
tematically studied and all qualitative results given below 
were supported by these systematics. Edge atoms, as we 
define them, are drawn in gray and are used to deter- 
mine edge-weighted vibrational density of states, where 
the edge-weight is ratio 



w 



dge 



i^edge i 



E 



(1) 



where (3-dimensional vector) is the component on atom 
i for eigenmode fj,. Hence w^^^g ~ 1 means complete lo- 
calization to edge, w^^ — J8t means "uniformly dis- 



tributed" mode. 



and w^^rfg. 



200 

would mean a mode lo- 



calized in the central region. Contributions of different 
directions are also analysed, analogously including weight 
of the components of the vibrational eigenvectors in 
either x-, y- or z-direction. 
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and similarly for y and z. Hence 
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edge 



(2) 



= 1 and 



wif measures how much 



of the mode is in the x-direction and ^^^^ measures how 
much of the mode at the edge is in the x-direction. 

Figure [5] shows the xx,xz and zz polarized Raman spec- 
tra of the ribbons, where e.g. xz stands for incident light 
polarization in x- and scattered light polarization in z- 
direction. The polarization pictures xx and zz show the 
same modes, but with different intensities. Polarization 
pictures that contain y-direction and out-of-plane modes 
yield tiny intensities and have been left out from figure [2l 
Orientation of the ribbons with respect to coordinate axis 
is shown in figure [TJ The most pronounced modes in figure 
[2] include breathing modes (intensive low-energy peaks 
in xx/zz spectra), G-band-related high-energy modes (in- 
tensive high-energy peaks in the xz and xx/zz spectra) 
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Fig. 3. Total vibrational density of states of zigzag, recon- 
structed zigzag (zigzag57) and armchair graphene nanorib- 
bons. Contribution of different directions in vibrational states 
are shown in different patterned areas. Vibrational density 
of states is broadened with normal distribution with a = 
45 cm~^, chosen to illuminate general changes in VDOS under 
reconstruction. 
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Fig. 4. As figure [3l but showing the edge weighted vibrational 
density of states of zigzag, reconstructed zigzag (zigzag57) and 
armchair graphene nanoribbons. 



and edge-localized modes (symbols (L) in xz and zz spec- 
tra). One of the main points is that only edge- localized 
modes undergo observable changes under reconstruction, 
others being within 2 cm^^ from corresponding zigzag rib- 
bon modes. These changes in edge-localized Raman active 
modes make the reconstruction visible. For zigzag57 the 
edge- localized mode of zigzag edge at 1550 cm~^ disap- 
pears in xz polarization spectrum in figure O and edge- 
localized triple-bond vibration becomes visible at energies 
around 2250 cm^^ in zz polarization spectrum of zigzag57 
ribbon, like in armchair ribbon. Intensity of the triple- 
bond vibration in reconstructed zigzag ribbon is smaller 
in our bond polarization theory than in armchair ribbon, 
mainly because of the differences in bond angles on the 
edge (the edge profile of zigzag57 is more linear). 

Couple of new Raman active modes can be noticed to 
appear for reconstructed zigzag edge, assinged with sym- 
bols (1) in XX and (2) in xz polarization spectrum in fi- 
gure [H Mode (1) at 1500 cm~^ includes stretching of edge 
bonds of pentagons. Mode (2) at 1600 cm^^ is mainly a 
longitudinal mode where first normal zigzag row of atoms 
in the zigzag57 ribbon vibrates as the edge-localized mode 
of zigzag ribbon. Neither of the modes (1) nor (2) is fully 
localized on the edge. 

Breathing modes remain unchanged, because they are 
collective vibrations which hold the edges internally un- 
changed during the vibration. Similarly, because G-band- 
related high-energy modes vanish towards the edge, re- 
construction has only minor effect on them. Raman active 
out-of-plane vibrations have much lower intensity than in- 
plane vibrations and are not therefore considered here. 

Next we compare more generally the vibrational pro- 
perties of the different ribbon edges. Figure [3] shows the 
total vibrational density of states, analysed with contribu- 
tions from different directions. The surprising and signi- 
ficant result is that zigzag57 ribbon's vibrational density 
of states resembles armchair ribbon's rather than zigzag 
ribbon's vibrational density of states. This means that for 
vibrational modes it is not the orientation of the under- 
lying honeycomb lattice that is important, but only the 
local structures of the edges. Comparison of changes bet- 
ween total vibrational density of states in figure [3] and 
edge- weighted vibrational density of states in figure 2] re- 
veals that the most significant changes are seen in edge- 
localized modes. Triple-bond vibrations in z-direction on 
the high energies appear as a separate peak around 2250 
cm~^ in figures [3] and [5] for both reconstructed zigzag rib- 
bon and armchair ribbon. The triple-bond edge modes 
in z-direction in the zigzag57 ribbon, somewhat deplete 
modes from the energy range 1700 — 2000 cm"-'^ of the 
G-band-related modes in z-direction. At the same time 
in the same energy range the number of edge-related in- 
plane vibrations in x-direction are increased due to recon- 
struction - induced edge stiffening. Similar effects are seen 
even in lower energies. Some edge-localized out-of-plane 
modes around 500 — 650 cm~^ are upshifted in energy 
for zigzag57 ribbon for the same reason. Even the energy 
of the lowest energy out-of-plane modes increase, which 
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can dramatically alter low-temperature properties (such 
as heat capacity) of the ribbons. 

Hence, due to edge stiffening, edge vibrations gene- 
rally up-shift in energy under reconstruction. This type 
of changes in mechanical properties can have an effect on 
applications such as sensors. The similarity of armchair 
and zigzag57 ribbon's VDOS is general observation, as 
the resemblance and also other qualitative features of the 
vibrational properties remained also for wider ribbons (up 
to = 40 A). 

To conclude, we have analysed the Raman spectra and 
vibrational density of states of the recently reported self- 
passivating edge reconstruction of zigzag ribbons. In the 
Raman spectra edge-localized triple-bond vibrations be- 
come visible like the corresponding mode in armchair rib- 
bon's spectra. These high-energy vibrations can be used 
to identify the reconstruction by concentrating the Raman 
measurement on the edge. 

The predicted intensity differences between triple-bond 
vibrations of the zigzag57 ribbon and the armchair ribbon 
could be used to separate the reconstructed edge struc- 
ture from the armchair edge structure with help of scan- 
ning tunneling microscopy images. Scanning tunneling mi- 
croscopy can also help in avoiding identification difhculties 
with ribbons that have combination of zigzag and arm- 
chair edges. For those ribbons the density of triple-bonds 
at the edges is lower and further intensity of the triple- 
bond vibrations is lower like in reconstructed zigzag rib- 
bon, which can lead to mix up between these structures. 

The rest of the visible Raman active modes of the rib- 
bons do not change for zigzag57 and the breathing mode 
has the same ribbon width dependence. Of course, in very 
narrow ribbons with majority of atoms on the edge, ener- 
gies or spatial nature of the G-band-related high-energy 
modes may still change. 
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